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In this paper we compare the relaxation in several versions of the Sznajd 
model (SM) with random sequential updating on the chain and square lattice. 
We start by reviewing briefly all proposed one dimensional versions of SM. 
Next, we compare the results obtained from Monte Carlo simulations with the 
mean field results obtained by Slanina and Lavicka . Finally, we investigate 
the relaxation on the square lattice and compare two generalizations of SM, 
one suggested by Stauffer and another by Galam. We show that there are no 
qualitative differences between these two approaches, although the relaxation 
within the Galam rule is faster than within the well known Stauffer rule. 

1 Introduction 

Applying the Ising model in the social sciences already has quite a long 
history jU Ej- In 2000 we proposed a simple model based on Ising spins 
[3], which was aimed at describing global social phenomena (sociology) by 
local social interactions (described by social psychology). The motivation to 
introduce the novel Ising spin dynamics was the so-called Social Validation - 
one of the most powerful phenomena that influence human decision j3] . The 
fundamental way that we decide what to do in a situation is to look at what 
others are doing. An isolated person does not convince others; a group of 
people sharing the same opinion influences the neighbors much easier. The 
crucial difference of our model, originally called USDF after the trade union 



1 



maxim " United we Stand, Divided we Fall" , compared to voter or Ising-type 
models is that information flows outward. The USDF model, later renamed 
by Dietrich Stauffer the Sznajd model (SM) has been modified and applied 
in marketing, finance and politics; for reviews see |H10IE]- 

Sociologically inspired models pose new challenges to statistical physics 
[13j. Therefore, simultaneously to applications, SM has been investigated 
from the theoretical point of view fOJEl- Moreover, numerous modifications 
of SM have been proposed [UJ HH [H] • In this paper we aim at gathering 
the models "under a common roof. Namely, we compare the relaxation in 
several versions of SM with random sequential updating on the chain and 
square lattice. 

2 One dimensional case 

The one- dimensional USDF model was defined in the following way [3]: 

1. In each time step a pair of spins Si and Si+i is chosen to change their 
nearest neighbors (nn), i.e. the spins and Si+2- 

2. If SiSi + i = 1 then Sj-i = Si and Si + 2 = Si - ferromagnetic rule. 

3. If SiSi + \ = —1 then Si-i = Si + \ and Si+2 — Si - antiferromagnetic 
rule. 

Two types of attractors are possible in such a model - ferromagnetic and 
antiferromagnetic steady states. Both are equally probable. 

Recently, attention has been paid to the second antiferromagnetic case. It 
was claimed that the antiferromagnetic rule could be considered to be quite 
unrealistic in a model trying to represent the behavior of a community. To 
avoid the unrealistic 50-50 alternating final state, new dynamic rules were 
proposed by Sanchez 0: 

1. In each time step a pair of spins Si and Si+\ is chosen to change their 
nearest neighbors (nn), i.e. the spins Si-i and Si+2- 

2. If SiSi+i = 1 then Si_i = Si and Si+2 = Si. 

3. If SiSi+i = -1 then S t -i = Si and S i+2 = S i+1 . 
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The model with these rules was studied by Monte Carlo simulations on a 
chain of length iV = 100. It was found that the antiferromagnetic final state 
was completely avoided and that the other two types of total agreement 
(ferromagnetic) final states could be achieved, with equal probability, from 
the same initial conditions. Moreover, the scaling exponent of the number 
of spins that changed their state with time was shown to be the same like in 
the original Sznajd model [TP] . 

This result is not very surprising, because both models (USDF |3] and 
Sanchez j^]) can be rewritten in a simpler form, namely: 

• USDF: Si_i takes direction of S i+ i and S i+2 takes direction of Si 

• Sanchez: Si-i takes direction of Si and Si + 2 takes direction of Si + \ 

In both dynamics a spin takes the direction of another spin for any value of 
SiSi+i- A similar idea was introduced by Ochrombel [TT] without qualita- 
tively affecting many of the results (except the presence of the phase transi- 
tions in two dimensions). In his approach a randomly chosen spin influences 
its neighbors, i.e. the neighbors get the same orientation. In one dimension 
this rule can be written as: 

• Ochrombel rule: Sj-i and Sj-fi take direction of Si. 

The observation that in all the described above dynamics only one spin 
influences other, makes the results obtained in ^2] more clear. In ^2] it was 
shown that one dimensional case USDF model could be completely reformu- 
lated in terms of a linear Voter model. 

Moreover, it is easy to notice that in the mean field approach all three 
models (USDF, Sanchez and Ochrombel) are equivalent. If we denote by N + 
the number of up spins and by N_ the number of down spins then we can 
write the magnetization defined as: 

N+-N- . 
m= — =p + -p-=p + -(l-p + )=2p+-l, (1) 

where p + is concentration of up-spins. 

In one time step, three events are possible: the magnetization increases 
by 2/N, decreases by 2/N or remains constant. For all three models (USDF, 
Sanchez and Ochrombel) these probabilities within the mean field approach 
for N oo are: 
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Number of analytical results such as mean relaxation time or distribution 
of waiting times have been obtained within such an approach [13J. 

As we have already said, in all three models only one spin influences its 
neighbor - the first neighbor for the Ochrombel and Sanchez modifications 
and the second neighbor for the USDF model. This is not very realistic from 
social point of view. Another possibility to avoid the antifferomagnetic final 
state, much more in the spirit of social validation, was briefly mentioned in 
[3| and investigated by the mean field approach in An analogous rule 
was proposed by Stauffer on the square lattice [Hj. From now on we call 
this rule Social Validation (SV): 

1. In each time step a pair of spins Si and Si + i is chosen to change their 
nearest neighbors (nn), i.e. the spins 6V i and S i+2 - 

2. If SiSi + i = 1 then = Si and Si + 2 = Si (social validation). 

3. If SiS i+ i = —1 then 6V1 = S^i and S i+2 = S i+2 (nothing happens). 

This model was studied by Slanina and Lavicka ^Hj using the mean field 
approach . They approximated the social network by the fully-connected 
network (the complete graph) of iV nodes. In such a case any two agents are 
neighbors and this is equivalent to the mean field approach. They were able 
to obtain results for the mean relaxation time: 



and deduce that the distribution of waiting times would have an exponential 
tail: 
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where r = t/N 2 = MCS/N is rescaled time, t is an elementary time step 
and MCS is one Monte Carlo step. 

All four presented above dynamics can be called "outflow" dynamics, 
since the information flows from the center spin (or spins) to the neighbor- 
hood. Of course we can investigate also "inflow dynamics" such as zero- 
temperature Glauber dynamics [To] : 

1. In each time step a spin Si is randomly chosen. 

2. If Si^Si+i = 1 then Si = S i+ i 

3. If Si^iS i+ i = —1 then Si = — Si with probability 1/2. 

Except for the original USDF model, all dynamics (Ochrombel, Sanchez, 
SV and Glauber) lead to full consensus on the chain with probability 1. A 
different situation is observed on the square lattice. A simple Ising ferromag- 
net has a large number of metastable states with respect to Glauber spin-flip 
dynamics fHl El ■ Therefore at zero temperature the system could get stuck 
forever in one of these states. There appears to be a nonzero probability 
that the square lattice system freezes into a stripe configuration. At T = 0, 
metastable states in this dynamics have an infinite lifetime that can prevent 
the equilibrium ground state from being reached. 

3 Monte Carlo results for one dimensional 
consensus rules 

In this section we compare the results obtained from Monte Carlo simulations 
on the chain of length N = 10 2 for all proposed one dimensional consensus 
models, i.e. Ochrombel, Sanchez, SV and Glauber dynamics. Next, we 
compare the results obtained from Monte Carlo simulations on the chain 
(N = 10 2 and N = 10 3 ) with the mean field results obtained in for SV 
dynamics. 

Probably the most natural way to investigate the relaxation process of the 
consensus dynamics is to look at the number of clusters in time. In all four 
dynamics the number of clusters monotonically decays as t" 1 / 2 . However, the 
differences between the dynamics can be observed if we look at the probability 
of reaching the final "all spins up" state (Fig.l) and the mean relaxation time 
(Fig. 2) as a function of the initial fraction of up-spins p+(0). 
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Figure 1: Probabilities of reaching the final "all spins up" state from a 
random initial state consisting of p+(0) up spins for Ochrombel (squares), 
Sanchez (stars), SV (triangles) and Glauber (circles) dynamics. 



It is seen that for three dynamics (Ochrombel, Sanchez and Glauber) 
the dependence between the probability P s t(all+) of the final "all spins up" 
state and the initial value of up-spins p+(0) is linear, while for SV the curve 
is S'-shaped (see Fig.l). This reminds of the phase transition obtained in the 
mean field approach J3| and for the simulations on the square lattice [Tlj . 

The differences between the dynamics are also clear if we look at the 
mean relaxation time < r st > as a function of the initial value of up-spins 
p + (0). For the Ochrombel, Sanchez and Glauber dynamics this dependence 
is well fitted by the quadratic function, while in the case of the SV model this 
dependence is bell-shaped (see Fig.2). Only quantitative differences can be 
observed between Ochrombel, Sanchez and Glauber dynamics - the relaxation 
under Glauber is much slower than under the other two. 

Let us now investigate the SV case more carefully. In Fig. 3 the proba- 
bilities of reaching the final "all spins up" state as a function of the initial 
value of up-spins p+(0) for two lattice sizes (L = 10 2 , 10 3 ) are presented. It 
can be observed that this dependence is not steeper for the larger lattice. 
Moreover, if we look at Fig. 4 we see that the mean relaxation time behaves 
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Figure 2: Mean relaxation times for an L = 100 chain from a random initial 
state consisting of p+(0) up spins for Ochrombel (squares), Sanchez (stars), 
SV (triangles) and Glauber (circles) dynamics. 
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Figure 3: Probabilities of reaching the final "all spins up" state from a ran- 
dom initial state consisting of p+(0) up spins for SV dynamics on the chain 
of length L = 10 2 and L = 10 3 . Clearly no phase transition is present on the 
chain. 



also the same for iV = 10 2 and iV = 10 3 . This suggests that there is no phase 
transition in one dimension, in contrast to the mean field results. 

However, the mean field results give the correct prediction on the distri- 
bution of the relaxation time. It was found that the distribution of waiting 
times had an exponential tail [T^]. Monte Carlo simulations confirm this 
prediction both on the complete graph ^3] and on the chain (see Fig. 5). 

4 Two dimensions 

Two dimensional models are much more realistic than one dimensional in 
describing social systems. Several possibilities of generalization to the square 
lattice were proposed by Stauffer et al. jHj. They presented the model on 
a square L x L lattice where again every spin can be up or down. Six 
different rules were introduced, but only two of them have been used in later 
publications: 
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Figure 4: Mean relaxation times (divided by a, to present results for L = 10 2 
and L = 10 3 on the same plot) from the initial state, containing randomly 
distributed p + = (m + l)/2 up-spins, to the ferromagnetic steady state. 
Obviously the mean field results are completely different from the simulation 
results. Moreover, comparison of the results for L = 10 2 and L = 10 3 clearly 
shows that no phase transition is present on the chain. 
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Figure 5: Probabilities of reaching the steady state in time larger than r on 
the chain for several values of the cinitial fraction of up-spins p. The distribu- 
tion of waiting times has an exponential tail with the exponent independent 
of p. 
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Figure 6: Dynamical rules of SM on the square lattice proposed by Stauffer 
(top) and Galam (bottom). In every time step a panel of four spins denoted 
by arrows is picked at random and influences spins denoted by x producing 
the configurations presented in the right-hand side pictures. 



• A 2 x 2 panel of four neighbors, if not all four center spins are parallel, 
leaves its eight neighbors unchanged (see Fig.l). 

• A neighboring pair persuades its six neighbors to follow the pair orien- 
tation if and only if the two pair spins are parallel. 

With both these rules complete consensus is always reached as a steady 
state. Moreover, a phase transition is observed - initial densities below 1/2 
of up-spins lead to all spins down and densities above 1/2 to all spins up 
for large enough systems [Hj. Galam (priv. comm. with Stauffer, described 
in J3J) showed that the updating rule of the one- dimensional SM can be 
transformed exactly into two dimensions in the following way (see Fig.l): 
the one-dimensional rule is applied to each of the four chains of four spins 
each, centered about two horizontal and two vertical pairs. 

Here we compare two rules in which a panel of four spins influences eight 
neighbors, i.e. Galam (Galam dynamics) and first of Stauffer's rules (Stauffer 
dynamics). It seems that Stauffer dynamics is more attractive from the social 
point of view j^]. Salomon Asch, in his famous experiment on conformity, 
found that one of the situational factors that influence conformity is the 
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Figure 7: Mean relaxation times from the initial state, containing randomly 
distributed p+(0) = (m + l)/2 up-spins, to the ferromagnetic steady state for 
Stauffer (left panel) and Galam (right panel) dynamics. Results are presented 
for the square lattice 30 x 30 (*) and 100 x 100 (o). The phase transition 
is clearly seen for both dynamics at p+{0) = 0.5. However, the relaxation 
under Glauber dynamics is faster than under Stauffer dynamics. 



size of the opposing majority. In a series of studies he varied the number 
of confederates who gave incorrect answers from 1 to 15. He found that 
the subjects conformed to a group of 4 as readily as they did to a larger 
group. However, the subjects conformed much less if they had an "ally". 
Apparently, it is difficult to be a minority of one but not so difficult to be 
part of a minority of two. 

On the other hand Galam dynamics is very simple for generalization on 
any regular lattice and was used to construct a two dimensional version of 
the so called two-component model [T8] . 

We have measured the probability of reaching the final "all spins up" 
state, as well as the mean relaxation time from a random initial state con- 
sisting of p+(0) up spins for Stauffer and Galam dynamics. We have found 
the phase transition for both dynamics - for p+(0) < 0.5 the "all spins up" 
state is never reached, while for p+(0) > 0.5 this state is obtained with prob- 
ability 1 (this result was obtained in [14). Moreover, critical slowing down 
is observed at p+(0) = 0.5 (see Fig. 7). For L — > oo we expect the 5(0.5) - 
function. 

We have also investigated the distribution of the relaxation time. In 
the mean field approach ^3] and on the chain (previous section) we have 
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found that the distribution of waiting times has an exponential tail with 
p-independent exponent (p = p + (0)). Results for the square lattice are pre- 
sented in Fig. 8. Under both dynamics the distribution of relaxation times 
has an exponential tail, but the exponent is p-dependent. Interestingly, the 
dependence between the exponent and the initial number of up-spins is iden- 
tical for both dynamics. 

Finally, we have measured the number of clusters N c in time under both 
dynamics. We have found that for the initial fraction of up-spins p = 0.5 
the number of clusters decays in time as a power law N c ~ t~ ' 5 (see Fig. 9), 
which is the same result as for the chain. Again, for any initial value of p the 
relaxation under Galam and Stauffer dynamics are quantitatively the same. 

5 Conclusions 

We have compared several one and two dimensional versions of SM. We 
have investigated the probability and the mean relaxation time of reaching 
the final "all spins up" state from the initial state, containing randomly 
distributed p+(0) up-spins. Moreover, we have measured the distribution of 
relaxation times and the decay of the number of clusters. We have found 
that on the chain Sanchez and Ochrombel modifications give qualitatively 
the same results as the zero-temperature Glauber dynamics, although time 
needed to reach the steady state under Glauber dynamics is much longer. In 
contrast the one dimensional dynamics proposed by Slanina and Lavicka in 
[T3*j . which we have called Social Validation dynamics is completely different. 
The difference is visible both in the probability and the mean relaxation 
time of reaching the final " all spins up" state from the initial state, containing 
randomly distributed p+(0) up-spins. For SV dynamics also analytical results 
were obtained ^3] and we were able to compare them with our simulation 
results. As often happens the results were remarkably different. In particular, 
no phase transition predicted by the mean field theory is observed in the one 
dimensional system. 

In the second part of this paper we have compared the results given by 
Stauffer and Galam generalizations of SM. This comparison seems to be quite 
important, since Stauffer generalization seems to be more attractive from 
social point of view, while Galam rule is much easier for generalization to 
other systems (in particular, it was used in the so called TC model [IE])- No 
qualitative difference has been found between these two dynamics, although 
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Figure 8: Probabilities of reaching the steady state in time larger than r on 
the square lattice for several values of the initial fraction of up-spins p for 
Stauffer (top) and Galam (bottom) generalizations of SM. The distribution 
of waiting times has an exponential tail with the exponent dependent on p. 
For both dynamics this dependence is identical. Critical slowing down is 
observed for p = 0.5 (right panels). 
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Figure 9: The number of clusters decays as a power law for both dynamics, 
time needed to reach the final state is shorter under Galam dynamics. 
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